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Abstract 

We develop a perturbative approach to the solution of the scalar wave equation for a large 

o : 

i ■ AdS black hole. In three dimensions, our method coincides with the known exact solution. We 
discuss the five- dimensional case in detail and apply our procedure to the Heun equation. We 
calculate the quasi-normal modes analytically and obtain good agreement with numerical results 

^— > , 

o for the low-lying frequencies. 
> 
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Quasi-normal modes describe small perturbations of a black hole at equilibrium [1-13]. A 
black hole is a thermodynamical system whose (Hawking) temperature and entropy are given 
in terms of its global characteristics (total mass, charge and angular momentum). In general, 
quasi-normal modes are obtained by solving a wave equation for small fluctuations subject to 
the conditions that the flux be ingoing at the horizon and outgoing at asymptotic infinity. One 
then obtains a discrete spectrum of complex frequencies. The imaginary part of these frequencies 
determines the decay time of the small fluctuations. On account of the AdS/CFT correspondence, 
the quasi-normal modes for an AdS black hole are expected to correspond to perturbations of the 
dual CFT. The establishment of such a correspondence is hindered by the difficulties in solving 
the wave equation. In three dimensions, the wave equation is a hypergeometric equation and can 
therefore be solved [13]. In five dimensions, the wave equation turns into a Heun equation which 
is unsolvable [15]. Numerical results have been obtained in four, five and seven dimensions [2, 14]. 

We discuss an analytic method of solving the wave equation in the case of a large black hole 
living in AdS space. The method is based on a perturbative expansion of the wave equation 
in the dimensionless parameter u/T H , where u is the frequency of the mode and T H is the 
(high) Hawking temperature of the black hole. Such an expansion is no trivial matter, for the 
dependence of the wavefunction on lv/Th changes as one moves from the asymptotic boundary 
of AdS space to the horizon of the black hole. The zeroth-order approximation is chosen to be an 
appropriate hypergeometric equation so that higher-order corrections are indeed of higher order 
in uj/T h . In three dimensions, this hypergeometric equation is exact. In five dimensions, we 
show that the low-lying quasi-normal modes obtained from the zeroth-order approximation are 
in good agreement with numerical results [2,14]. The first-order correction is also calculated. 
Our method is generalizable to higher dimensions where the number of singularities of the wave 
equation is larger. 

The metric of a <i-dimensional AdS black hole may be written as 

ds 2 =(^- 2 +l- U -^f] + 7 ^ . + r^ d 2 (1) 

where R is the AdS radius and M is the mass of the black hole. For a large black hole, the metric 
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simplifies to 

= ( 

R 2 r .d-3 
The Hawking temperature is 



( r 2 _ <£d-\M \ 



r »=— # < 3 » 



where r fc is the radius of the horizon, 



The scalar wave equation is 



R d-3 



i/(d-i) 

(4) 



-^=d A y/gg AB d B ^ = m 2 $ (5) 

We are interested in solving this equation in the massless case (m = 0) for a wave which is 
ingoing at the horizon and vanishes at infinity. These boundary conditions yield a discrete set 
of complex frequencies (quasi- normal modes). 

We start with a review of the three-dimensional case where the assumption of a large black 
hole is redundant and the wave equation may be solved exactly. Indeed, in three dimensions 
(d — 3), the metric reads 

1 R 2 dr 2 

ds2 = t^ (r 2 - r D df2 + t~9 — 2Y + r2dx2 ( 6 ) 

R 2 v ' (r 2 -rl) 

independently of the size of the black hole. The wave equation is 

W~r dr i 1 ~ 7 2 ) ^) - r^rf^ + V 2 ^ ~ ^ = ° (?) 
The solution may be written as 

r 2 

$ = e i ( wt -P a: )*(y), y = -| (8) 

where * satisfies 

y 2 (y - 1) ((y - 1)*')' + ^ 2 + f y(y - 1) + I^ 2 (y - 1)* = (9) 
in the interval < y < 1, and we have introduced the dimensionless variables 



Two independent solutions are obtained by examining the behavior near the horizon (y — > 1), 

* ± „(l_ y )±tt (11) 

where ^ + is outgoing and is ingoing. A different set of linearly independent solutions is 
obtained by studying the behavior at large r (y — > 0). We obtain 



*~y h± , h± = -±-VT+W (12) 

In the massless case (m = 0), we have h + — 1 and /i_ = 0, so one of the solutions contains 
logarithms. For quasi-normal modes, we are interested in the analytic solution. This may be 
written in terms of a Hyp ergeo metric function as 

*(y) = y(l - yfiF^l + i(u+p),l + i(u - p); 2; y) (13) 

Near the horizon, this solution is a mixture of ingoing and outgoing waves. Using standard 
identities involving Hypergeometric functions, we obtain 

V ~ A(l - y)-** + B(l - y)** (14) 

where 

A = T jm B = n-^) (15) 

r(l +i(u + p))T(l +i(u-p)) ' T(l-i(u+p))r(l-i(u-p)) v ; 

as y — > 1. Since ^ must be a linear combination of and we deduce 

^ = A^_ + J B^+ (16) 

on account of eq. (11). For quasi- normal modes, we demand that \& be purely ingoing at the 
horizon, so we need to set 

B = (17) 
with B given in (15). The solutions to this equation are the quasi-normal frequencies. Explicitly, 

Cj = ±p — in , ra=l,2,... (18) 

a discrete set of complex frequencies with negative imaginary part, as expected [2]. Notice that 
we obtained two sets of frequencies, with opposite real parts. 



In five dimensions (d — 5), the wave equation (5) reads 



ld r r 5 (l-d\d r Q ^ 4x a 2 $ - % V 2 $ - m 2 i? 2 $ = (19) 



Let 



and change parameter r to y, 



The wave equation then becomes 



$ = e^*-^)^(r) (20) 



( 21 ) 



where 



y 3 (i - y 2 ) Q(i - y 2 )*') + y y* + f y(i - y 2 )* - X(i - y 2 )* = o (22) 

Near the horizon we obtain in- and out-going waves 

#± ~ (l - y) ±li/4 (24) 



v(y)=y 2 (i-yr* j/A (- 



At large r (y — > 1) there is only one analytic solution for m = (the other one contains 
logarithms). It may be written as 

)-u>/4 
F{V) (25) 

where we isolated the singularities at y — 0, ±1. We divided (1 + y) by 2 in order not to obtain 
a contribution from this factor at the horizon (y — > 1). This introduces an irrelevant constant, 
but one ought to be careful with Co dependent constants in perturbation theory (expansion in 
Co). It should also be pointed out that we could have chosen the exponent +Co/A instead of —Co /A 
at the y — — 1 singularity (eq. (25)). This does not alter the quasi-normal modes, but in the 
perturbative expansion, each choice only produces half the modes. The two sets of modes have 
the same imaginary parts, but opposite real parts, as we shall see. 
F is the Heun function satisfying the equation [15] 

F „ + f 3 + l-jfl/2 + F , + (2 - (1 + ,>/4ft, - , p = Q 

\y y- 1 y + 1 / y(y - 1) 



where 



4 4 4 



To find the behavior of F near the horizon, let us write the Heun equation as 

(Ho + Jii)F = (28) 

where 



H = x(l-x)^ + (2-l T i^-(3-^^)A_l ((2 -l±i^_ g ) 



(29) 



and we changed variables to x = y 2 . We shall treat Hi as a perturbation. It has been chosen so 
that the perturbation series is an expansion in Cj ~ uj/T h . Expanding the wavefunction, 

F — F + Fi + ... (30) 

the zeroth-order equation 

n F = (31) 
has analytic solution the Hypergeometric function 

F (x) = 2 F 1 (1 + (-l±i co + Vg)/2, 1 + (-±±* £ - Jq)/2; 2-^u;x) (32) 

Its behavior at the horizon (x — > 1) is 

Fo(rr) ~ Ao + Bo(l-x)^ 2 (33) 

where 

r(2- ±fi£)r(iw/2) 



A ° r(i - (±=* £ + ^)/2)r(l - (±=* * - V?)/2) 
r(2- ifiw)r(-iw/2) 
r(l - (±±^ + v /g)/2)r(l-(±±^ - ^)/2) 



Bo r(i - (±±* & + ^)/2)r(i - (ip u> - /f " (34) 

An approximation to the quasi-normal modes is obtained by setting 

B = (35) 
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The solutions to this equation are 

i^u±y/q = 2n, n = l,2,... (36) 

For p — 0, they are shown on Table 1, column (a). They all have negative real parts. The set of 
modes with positive real parts may be obtained by replacing the factor (lia)-^/ 4 with (I±l)+ W / 4 
in eq. (25) (capturing the alternative behavior of the wavefunction near the singularity y — — 1). 
The low- lying frequencies are close to the exact values obtained by numerical methods (Table 1, 
column (c) [14]). For first-order perturbation theory, we need the expansion of B in a), 

To calculate the first-order correction, we also need another linearly independent solution of the 
zeroth-order wave equation. We shall use 

r w (x>)dx> 

Go(x) = F (x)J (38) 

where W is the Wronskian 

Wo = F G' - G Fq = x~ 2+(1 -^/\l - X y l+i "/ 2 (39) 
The equation satisfied by the first-order contribution to the wavefunction F is 

H F 1 = -HiF (40) 
The solution (satisfying Fi(0) = 0) may be written as 

Fl(x) = Go(x) [' f.wh,FoW^ _ r a^yn^w (41) 

Expanding in a), we obtain from eqs. (32) and (38), respectively, 

W = L J +... , G (x) = — + ... (42) 

iujx/2 x 

The first term in F x (eq. (41)) does not contribute at lowest order. After some algebra, we arrive 
at 

F 1 (x) = C±fcuF (x) + ... (43) 
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where 

7T 2 

= l + ln2 (44) 



Notice that the divergences from the contributing terms in the integral in eq. (44) at both ends 
(x = 0, 1) all cancel each other, resulting in a finite expression for C, as required for the validity 
of our perturbative expansion. To first order, the behavior at the horizon is 

F(x) « F (x) + F 1 (x) ^A, + B 1 (l - xf t2 (45) 

Using eqs. (37) and (44), we obtain 

B 1 = B + 1 -^C =— l — {l + ln2±f £ + ...} (46) 

At this order, we obtain an approximation to the lowest quasi-normal frequency by setting B\ = 
in eq. (46). We find 

^i = -^^ = -2.89(1 + (47) 

which is a good approximation to the exact result CJ\ = ±3.12 — 2.75i obtained by numerical 
techniques [2,14]. (The mode with positive real part is obtained by replacing the exponent 
— Co /A by +lu/4 at the singularity y — — 1 (eq. (25)), as already explained). Higher-order modes 
may be obtained by calculating higher-order perturbative corrections. At the nth perturbative 
order we obtain an nth order polynomial in uj whose roots are an approximation to the n lowest 
quasi-normal frequencies. 

Our method may be straightforwardly extended to higher dimensions even though the wave 
equation possesses more (complex) singularities than the Heun equation in five dimensions. It 
would be interesting to investigate the implications of our perturbative approach to the AdS/CFT 
correspondence. 
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n 


(a) 


(b) (c) 




1 


±1.19 - 3.03i ±2i 


*9 - 2.89i ±3.12 - 


2.75i 


2 


±5.01 - 7.65i 


±5.17- 


A.76i 


3 


±9.01 - 11.78i 


±7.19- 


Q.77i 


4 


±13.01 - 15.84i 


±9.20 - 


8.77i 


5 


±17.00 - 19.87i 


±11.20- 


10.77i 


6 


±21.00 - 23.90i 


±13.21 - 


12.78i 


7 


±25.00 - 27.92i 


±15.21 - 


14.78i 


8 


±29.00 - 31. 93i 


±17.21 - 


16.78i 


9 


±33.00 - 35.93i 


±19.21 - 


18.78i 


10 


±37.00 - 39.94i 


±21.21 - 


20.78i 



Table 1: Quasi-normal frequencies in d — 5: ^aj zeroth-order untruncated approximation 
(eq. (35)), f&j first-order approximation (eq. (47)), (c) numerical results [14]. 
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